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The non-collinear spin-spiral density wave of the uniform electron gas is studied in the framework 
of Reduced-Density-Matrix-Functional Theory. For the Hartree-Fock approximation, which can be 
obtained as a limiting case of Reduced-Density-Matrix- Functional Theory, Overhauser showed a long 
time ago that the paramagnetic state of the electron gas is unstable with respect to the formation of 
charge or spin density waves. Here we not only present a detailed numerical investigation of the spin- 
spiral density wave in the Hartree-Fock approximation but also investigate the effects of correlations 
on the spin-spiral density wave instability by means of a recently proposed density-matrix functional. 
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I. INTRODUCTION 



For many decades, the uniform electron gas has 
served as the model for the description of many-particle 
systems^. However, the determination of its ground 
state, without any symmetry assumptions, still remains 
a challenge. Specific symmetries for the fully correlated 
uniform electron gas have been investigated using Monte 
Carlo methods 2 -! 3 -. These studies focus mostly on broken 
spatial symmetry, i.e. , Wigner crystallization, or broken 
global spin symmetry. 

For the electron gas with constant electron density and 
uniform spin-polarization, the ground-state energy is an- 
alytically accessible in the Hartree-Fock approximation. 
Overhauser showed in his seminal work 4 ^ that within 
the Hartree-Fock approximation the aforementioned ho- 
mogeneous ground state exhibits an instability w.r.t. the 
formation of charge and spin density waves. Wigner 
crystallization within Hartree-Fock has been investigated 
in Ref. |g. Only recently the combined local spatial- 
and spin-symmetry breaking of the Hartree-Fock ground 
state has been studied using a Monte Carlo method which 
optimizes the ground-state energy in the space of single 
Slater-determinants 7 . However, this study still remains 
in the regime of collinear spin polarization. 

In the present work we investigate the case of lo- 
cal spin symmetry breaking, specifically a non-collinear 
spin-spiral symmetry. We employ Reduccd-Density- 
Matrix- Functional Theory both in the limiting case of the 
Hartree-Fock approximation as well as for the correlated 
electron gas using the recently proposed density-matrix- 
power functional^. 



II. THEORETICAL FRAMEWORK 

A. Reduced-Density-Matrix-Functional Theory 

The basic variable in Reduced-Density-Matrix- 
Functional Theory (RDMFT) is the one-body-reduced 
density matrix (1-RDM) defined by 

1<ya , (r;r') = Ti N [b^l (/) & (r) } (1) 

where D is the zero-temperature statistical operator of 
an ensemble of iV-electron states 

d^Y,^M)(^\ with E^ 2 = 1 > ( 2 ) 

i i 

where i/ , i( r ) and ipo-(v) are fermionic creation and anni- 
hilation operators, respectively. The 1-RDM is a Hermi- 
tian operator in the single-particle Hilbert space and can 
be represented by its spectral decomposition 

7 (r;r') =^ni*i(r)$t(r'), (3) 

i 

where the eigenvalues nt are called occupation numbers 
(ON) and the corresponding single-particle Pauli-spinor 
eigenstates $,(r) = (ip^(r) ,(pii(r)) are referred to as 
natural orbitals (NO). It was shown by Gilbert that 
the iV-particle ground state is a unique functional of the 
ground state 1-RDM, i.e. , |*^) = |*^[7 gs ]) . Therefore 
the ground state energy for a system of N interacting 
electrons moving in an arbitrary but fixed (possibly non- 
local) external potential V is also a functional of the 1- 
RDM: 

£v[7 gs ] = «[f s ]|Wy|<[f s ]) , (4) 
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where Hv — T + V + W + Ei on is a generic interacting 
many-body Hamiltonian with kinetic energy T, external 
potential V, electron-electron interaction W, and a con- 
stant energy contribution -Eion from the degrees of free- 
dom that are not treated quantum mechanically. 

The ground-state-energy functional can be decom- 
posed into the following components 



Ey[y] = T[ 7 ] + V[j] + W[y] + E ion , 



(5) 



with the kinetic energy (atomic units are used through- 
out the paper, and the superscript "gs" is omitted for 
brevity) 

T M =E ! im Vv 7 ..(r;r'), (6) 

and the energy contribution due to the external potential 



V[~l]=Y J Jd 3 rV(r) lr7rT (r;r) 



(7) 



Here we are assuming a local, spin independent exter- 
nal potential. The Hohenberg-Kohn theorem of Density- 
Functional Theory (DFT) proves a one-to-one mapping 
between the ground-state density and the N-particle 
ground state, considering only local external potentials. 
However, in RDMFT the Gilbert theorem ensures a 
one-to-one correspondence between the ground-state 1- 
RDM and the N-particle ground state by considering the 
broader class of non-local external potentials. This also 
implies the one-to-one mapping between a local poten- 
tial and the ground-state 1-RDM. Note that in contrast 
to usual Kohn-Sham DFT all single-particle contribu- 
tions to the ground state energy Ey are explicitly given 
in terms of the ground state 1-RDM. However, the inter- 
action energy 



W[j) 
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(8) 



is only known explicitly in terms of the ground-state pair 
density 

(^Wl^Cn)^^)^^)^^)!^^]) . (9) 

The basic idea of RDMFT is to extend the domain 
of the ground-state-energy functional in Eq. to all 
ensemble- iV-representable 1-RDMs [as defined in Eq. |T])] 
and then employ the variational principle in order to find 
the ground-state 1-RDM as well as the ground-state en- 
ergy corresponding to a fixed external potential V. The 
necessary and sufficient conditions for a 1-RDM to be 
ensemble- iV-representable areii: 

E rii = N and < m < 1, (10a) 

i 

E|d 3 r^ CT (r)^(r)=%. (10b) 



In order to apply RDMFT in practice we need to approxi- 
mate the functional dependence of the pair density on the 
1-RDM. Since we want to study the spin-spiral density 
wave (SSDW) instability in the uncorrelated (Hartree- 
Fock, HF) and the correlated regime, we focus on the 
so-called density-matrix-power functional introduced in 
Ref. § 

P^a 2 [H(ri, r 2 ) =2^1*1 {rv, n) 1<? 2 * 2 (r 2 ; r 2 ) 

for 0.5 < a < 1. Here the power of the 1-RDM has to be 
read in the operator sense, i.e. 



7°(r; rO = £"?*<to*i CO- 



(12) 



As limiting cases it contains both the uncorrelated HF 
approximation (for a — 1) as well as the correlated 
Miillcr or Buijsc-Baercnds functional (for a = 0.5) 12 ' 13 . 
Also, it was recently shown^ that the power functional 
yields good correlation energies for the unpolarized uni- 
form electron gas. 



B. The Overhauser Instability of the uniform 
electron gas 

The system under investigation is the uniform electron 
gas (UEG) in three dimensions, i.e. , a gas of interacting 
electrons subject to an external potential induced by a 
uniformly distributed positive background charge. Over- 
hauser has proved that the true HF ground state does 
not correspond to a homogeneous electron density (al- 
though there are solutions to the HF equations where 
the symmetry is not broken), since the HF energy can 
be lowered by forming a charge density wave (CDW) or 
spin density wave (SDW) 5 . As an explicit example he 
assumed, in addition to the regular HF potential V a ki a 
potential gu in the HF Hamiltonian that couples plane 
waves of opposite spin whose momenta differ by q: 



n HF = 



ha 



Overhauser demonstrated that with the ansatz 




sm 



cos 



q r 



q r 



(14a) 
(14b) 



the HF self-consistent equations are transformed into a 
set of equations relating the orbital angles 9k, the poten- 
tial gk and the regular HF potential V a k- Note that the 
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generic single-particle index i here has been replaced by 
the joint index i — >{(& = 1, 2) , k}. After taking the ther- 
modynamic limit ( i.e. , the volume and the number of 
particles N are taken to be infinity such that |y remains 
constant), these equations read (cf. Ref. 1 141) 

'd 3 fc' 4tt 



Tfe-f 



(2tt 
x < cos 



k - k'\ 
2 I ^h' 



n 2 k' sm 



(15a) 



v. 



ifc+f 



d 3 fc' 
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x < n^.' sm 



fe - ft/ 

2 / ^fe' 



n 2 fc' cos 



2fffc 



d 3 fc' 4tt 
(^L^f^' 



(15b) 

n 2 fe'}sin(0 fe /). (15c) 



The r.h.s. of Eqs. (| 1 5(1 implicitly depends on q via the 
ribfc and the 9k- The ribfc are the occupation numbers 
(either or 1) of the orbitals <!>{,/- which comprise the HF 
ground-state Slater-determinant and specify the Fermi 
surface (the boundaries of the integration) in Eqs. (|15a|) 
- (|15cp . The orbital angles 6k on the other hand are given 

by 



tan(0 fe ) =- 



(fc-|g) 
2 



U T , 



-V, 



J.fc+4 



(16a) 

(16b) 
(16c) 



Note that the origin in momentum space is shifted by q/2 
compared to the definitions in Ref.ll4l The energy contri- 
bution due to the pairing potential favors a hybridiza- 
tion of spin-up and spin-down plane waves differing by 
q in their momenta. The orbital angles 9k introduced 
in Overhauser's ansatz Eq. (fT*!)) describe this hybridiza- 
tion. Another way of looking at the orbital angles Ok is 
to consider them, together with the angles </>(r) = q-r, 
as angles defining a rotation in spin space represented by 



U(r;k) 



cos(ifl fc ) e-** r - sin(i0 fe ) e~^ r 



(17) 



where a y l z are Pauli matrices. The orbitals of 
Eq. (|14ap [ (|14b[) ] can then be thought of as being con- 
structed by transforming pure spin-up [spin-down] plane 
waves in spin space according to the rotation Eq. (fTT)) . 
First the plane wave is rotated around the y-axis by 




FIG. 1: (Color online) The effects of the spin rotation U(r ; k) 
on pure spin-up (dashed arrow) or pure spin-down (solid ar- 
row) natural orbitals (plane waves) for two momenta fci/2- 
The angle 9k specifies the cone on which the spin is rotating. 
The position on the cone is given by the angle (j>( r ) = Q r > 
which is the same for all natural orbitals. 



an angle 9k, i.e. , an angle depending on its momen- 
tum. Then it is rotated around the z-axis by an angle 
<f)(r) = q r which is the same for all plane waves, inde- 
pendent of the wave vector, but depends on the spatial 
position (see Fig. []}. With this consideration it is clear 
that the angle Ok has to be restricted to the interval [0, ir] 
in order to assign a unique azimuthal rotation angle. 

In previous studies within RMDFT-^ it was assumed 
that the 1-RDM exhibits the symmetries present in the 
Hamiltonian, i.e. , the NOs are pure spin-up(down) plane 
waves, while here we use orbitals of the form of Eq. (jT5J) 
as NOs for our RDMFT treatment of the UEG. The spin- 
spiral wave vector q and the angle 9k will be treated as 
variational parameters for the NOs. It can easily be ver- 
ified that the NOs of Eq. (fT^r) form a complete and or- 
thonormal set and that the corresponding electron den- 
sity p = 3/(47rr 3 ), given in terms of the Wigner-Seitz ra- 
dius r s , is still spatially constant. The magnetization of 
the UEG is defined by 



m(r) =- £X;<tf|#/(r)ow&(r)|¥) 



(18) 



Kj n {r;r) 
X7 Ti (r;r) 

i{7n( r ; r ) -7u0;r)}, 



and varies in space as 

'A cos(q-r) N 
m(r) = — | Asin(q-r) 
B 



A 
B 




{nik - n 2 k}sm(9 k ) ■ 



{nik - n 2 k}cos(0 k ) 



(19a) 

(19b) 
(19c) 



i.e. , the x- and y-components of the magnetization ro- 
tate in space along the direction of q with a periodicity 
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given by the wavelength q — \q\. This geometry of the 
magnetization is usually referred to as SSDW 28 . 



III. NUMERICAL IMPLEMENTATION 

Having chosen a functional and having made an ansatz 
for the NOs, we minimize the functional for the ground- 
state energy. The functional depends on n b k , 8k and the 
spin-spiral wave vector q. The contribution Ei on coming 
from the uniform positive background charge cancels ex- 
actly the classical contribution of the interaction energy, 
since the density is constant. Accordingly the energy per 
electron reads 

e a [n b ,0](q) = t{n b ,9](q) - w al {n b ,8] - w a2 {n b ,8] , (20) 

with the kinetic energy per electron 

{2-Kf 

q-k(n lk - n 2fc )cos(6> fc )} + y , (21) 



1 f d 3 k r 
t[n b ,8](q) =— J— ^{{ni k + n 2k )> 



the energy contribution from exchange-like terms of or- 
bitals with the same b (intra-band exchange) 



w al [n b ,( 



1 ffd 3 k 1 d 3 k 2 4tt 
Yp 



[2nf (fe!-fc 2 ) 2 

x{(nik x n lk2 ) a + {n 2kl n 2 k 2 ) a } 

■> I 8ki — 8k 2 
X cos 1 



(22) 



and the energy contribution from exchange-like terms of 
orbitals with opposite b (inter-band exchange) 



w a 2[n b ,8] 



1 /7d 3 fcid 3 fc 2 



-in 



{2nf (fc!-fe 2 ) 2 



2p 

X {(«ifc 1 «2fc 2 ) Q +(n 2 k 1 nik 2 ) a } 

> / 8k x — 8k 2 
x sin ' 



(23) 



We assume that the symmetry is only broken along the 
direction of q which is chosen to be parallel to the z- 
axis. Accordingly we can use cylindrical coordinates in 
momentum space, i.e. , n bk = n bkpkz and 8 k = 8 kpkz . We 
also use the following additional symmetry assumptions 



nik > n 2k , 



n bkp - kz =n bkpkz 
8k p ±\k z \ T a kp \ kz \) , 



(24a) 
(24b) 



with < Ofe < 1. In this way we guarantee that the en- 
ergy gain in the part of the energy which explicitly de- 
pends on q is maximized. The z-component of the mag- 
netization vanishes under these symmetry assumptions 
(planar spiral). 



The configurations 



^ =9(|fe - e z k f \ - kf) + Q(\k + e z k f \ - k f ) . 



,PM 



>2k =0 , a£ M = 1,Q= 2fc/e 2 
and 



FM 
l lk 



=e( k ~ 2 1/3 k f 



nL M -0, a™ = 0, 9 = 0, 



(25a) 



(25b) 



which are compatible with Eq. (|24[) . correspond to the 
non-magnetic (usually in this context called paramag- 
netic [PM]) and ferromagnetic (FM) state of the UEG 
within HF, respectively. 

When discretizing the integrals of Eqs. (f2"Tj) - (|23[) we 
assume that the ONs n b k and the angles ak are constant 
within annular regions in fc-space 



{fe | k;~ <k p < k p + ; k\- <k z < K+ }. (26) 



Then the discretized energy contributions are 

q 2 

t[nu,0i](q) = ^ n b tBKli + — 

bi 

- g^(nii - n 2i ) 003(8,) DQI t , (27a) 

i 
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where the integral weights are given by 



DKL- = 



DQI, 



1 



8tt 2 p 



£2, 



: i i d.kpdk z (kpk z ) 



(27b) 



(27c) 



(28a) 



(28b) 



nvi 1 fff dk pl dk zl d<j>i fff dk p2 dk z2 d(j> 2 
DXI " ' ; ~ ' 1 1 ^~JJJ Qi (2vf 



4irk pl k p2 



k li + k p2 +( k zi ~ k z2 ) - 2k pl k p2 cos(0i - 02 ) 



(28c) 



The integrals (|28a[) and (|28b[) are readily solved and the 
integrals (|28c|) can ultimately be reduced to elliptic in- 
tegrals, which are numerically accessible with high ac- 
curacy. Since the momenta are treated as continuous 
variables we stay in the thermodynamic limit. Thus all 
energies obtained numerically are variational. The error 
introduced by the discretization is solely due to the as- 
sumption that the n b k and 8k are constant within the 
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FIG. 2: (Color online) Hartree-Fock total energy per electron 
of the SSDW state as function of the spin-spiral wave vector 
q at various r„. The value e(g = 2k f) is subtracted in order 
to emphasize the behavior of the minimum at different den- 
sities. For increasing density (decreasing r s ) the minimum 
shifts to higher values of q opt and the energy gained against 
the paramagnetic state by forming a spin spiral decreases. 




-0.047 



-0.048 



FIG. 3: (Color online) Dependence of the energy per electron 
on r 3 , within the Hartree-Fock approximation, for the para- 
magnetic, ferromagnetic and SSDW phase in the region of 
the paramagnetic-ferromagnetic crossover. The inset shows 
the behavior of the amplitude A, defined in Eq. (|19b|l . at the 
optimal spin-spiral wave vector. 



elementary volume elements tti and can systematically 
be reduced by increasing the number of discretization 
points. 

After having discretized the problem, the minimiza- 
tion of the energy functional of Eq. (j2"0|) becomes a high- 
dimensional optimization problem. We use a steepest 
descent algorithm for the minimization and ensure that 
the constraints, Eq. (TIT)]) , are satisfied during the mini- 
mization process. Starting from some initial 1-RDM and 
some initial discretization in momentum space the energy 
is minimized for a fixed spin-spiral wave vector q. Then 
the discretization is refined in those regions of momen- 
tum space where the and/or the cii show the largest 
variations. The minimization on the refined momentum 
space mesh starts from a re-initialized 1-RDM in order 
to prevent dependencies on the result of the minimiza- 
tion on the coarser grid. Finally we compare the total 
energies at different q in order to determine the optimal 
spin- spiral wave vector q t for various densities. 



IV. RESULTS 

A. Hartree-Fock 

We first use our numerical implementation to investi- 
gate Overhauser's SSDW state in the HF approximation, 
i.e. , the density-matrix-power functional with a = 1. 
From the considerations in Eqs. (|25[) we see that it is 
sufficient to minimize w.r.t. a 1-RDM whose ONs are 
only non-zero for orbitals with b—l and \q\ € [0, 2kf] 
since both the paramagnetic and the ferromagnetic HF 
solutions are accessible under these conditions. The min- 
imization at q — and q — 2k f yields exactly the ONs nu 
and angle parameters given in Eqs. (|25b[) and (|25a|l . 



respectively. Therefore we can read the total energy per 
particle as function of the spin-spiral wave vector in the 
following way: e(q — 0) is the energy of the ferromagnetic 
state, e(q = 2k f) corresponds to the energy of the para- 
magnetic state. For intermediate values, < q < 2k j, 
e(q) corresponds to a SSDW configuration with m z = 
(planar spiral). Overhauser's statement can then be ex- 
pressed as d q e(q)\ q _ 2k > 0, i.e. , the paramagnetic con- 
figuration is unstable w.r.t. the formation of a SSDW. 

In Fig. [5] we show the dependence of the total energy 
per particle on the spin-spiral wave vector q for vari- 
ous densities. Consistent with Overhauser's proof, the 
derivative of e(q) is positive at q = 2kf. It is clear from 
Fig. [2] that the optimal spin-spiral wave vector moves 
away from the paramagnetic configuration (q — 2kf) as 
the density decreases. Furthermore the difference be- 
tween the total energy at the minimum and the total 
energy at q — 2k f increases with increasing r s , i.e. , the 
instability is more pronounced at lower densities. Below 
some critical density, however, the ferromagnetic state 
(q = 0) becomes the most stable solution. This is not 
in contradiction with Overhauser's statement since the 
spin-spiral state is still lower in energy than the param- 
agnetic state. A comparison of the energy per electron in 
the paramagnetic, ferromagnetic and SSDW phase is de- 
picted in Fig. [3l We provide results for the non-collincar 
magnetic states of the UEG in order to extend the pic- 
ture given in Ref. 0. It seems that the gain in energy 
by forming a collinear SDW/CDW state as presented in 
Ref. lD is larger compared to the energy gain by forming a 
SSDW. This is consistent with the qualitative argument 
already given by Overhauser, that the superposition of a 
left- and right-rotating SSDW yielding a collinear SDW 
will increase the gain in energy*. 

To describe the resulting behavior of q op t(r s ) we pro- 
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FIG. 4: (Color online) Dependence of the Hartree-Fock opti- 
mal spin-spiral wave vector q op t on the density, given by r s . 
The proposed approximation by a simple scaling law Eq. (|29p 
is shown as the dashed line. 
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FIG. 5: (Color online) Hartree-Fock total energy per electron 
as function of the spin-spiral wave vector q at the density 
corresponding to r s = 5.0. The data sets represent results at 
different discretizations. The dashed horizontal line visualizes 
the analytic value of the paramagnetic ground state energy. 
The optimal spin-spiral wave vector is q opt ~ l.Qkf. 



pose a simple, empirical scaling law for the optimal spin- 
spiral wave vector 



q op t(r s ) = 2kf[ 1 - ( — 



(29) 



where r$ w 5.7 and /3 pa 0.2. The proposed scaling behav- 
ior of q opt reproduces the numerical data very accurately 
as can be seen in Fig. [4] It should be emphasized that we 
do not find any optimal spin-spiral wave vector q opt < kf. 
Note that for densities close to the transition to the fer- 
romagnetic state the optimum wave vector q opt can be 
quite different from 2k f while for higher densities it is 
very close to this value. 

The effect of the refinement of the discretization in mo- 
mentum space is shown in Fig. [5] By sampling and 



FIG. 6: (Color online) Hartree-Fock single-particle dispersion 
(k p — 0) at r s — 5.0 for the paramagnetic (q > 2kf), ferro- 
magnetic (q = Okf) and the SSDW state (q op t = 1.6fc/). The 
single-particle energies are shifted such that the dashed hor- 
izontal line corresponds to the Fermi energy for all q. The 
difference between the two symmetric minima corresponds to 
the spin-spiral wave vector q. The paramagnetic dispersion 
may also be viewed as a spin-spiral dispersion with the ori- 
gin in momentum space shifted by ±q for the different spin 
channels [cf. Eqs. (55} ]. 



more often in regions of higher variations we both lower 
the energy and reduce the numerical noise in e(q). The 
convergence of the total energy can be inferred from the 
values e(q = 2fc/) at different discretizations and compar- 
ing to the analytic paramagnetic energy. For the case of 
r s — 5.0 we obtain a spin-spiral energy that is lower than 
the analytic paramagnetic energy at the optimal value of 
the spin-spiral wave vector. At higher densities (lower 
r s ) the energy gain by forming a SSDW is lower, so we 
would need a very fine discretization to obtain numeri- 
cal results lower than the analytic paramagnetic energy. 
However, considering the numerical value of the param- 
agnetic energy at the same discretization is sufficient to 
demonstrate the instability w.r.t. a SSDW formation be- 
cause the computed energies are variational as discussed 
in Sec. IIIIl In order to determine the dependence of the 
optimal spin-spiral wave vector q opt on the density, we 
therefore refine the momentum space discretization until 
q opt is converged. 

For our numerical results we have verified that the ONs 
and the angular parameters satisfy Overhauser's self- 
consistent equations ([To]) and (flU)) by iterating them only 
once. The difference between the angles in the occu- 
pied regions before and after the iteration is numerically 
zero for all values of q. This means that choosing a spin- 
spiral wave vector we can always find a solution of the 
self-consistent equations derived by Overhauser. Since 
the total energy does not depend on the in regions 
where nu = 0, one self-consistency loop furthermore fixes 
the angles in unoccupied regions of fc-space because 
they appear only on the left-hand side of Eqs. (|16p . 
This is necessary to construct the proper HF dispersions 
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FIG. 7: (Color online) Orbital angles 9(k p = 0,k z ) for vari- 
ous densities, specified by r s , at the optimal spin-spiral wave 
vector. The horizontal dashed line corresponds to the orbital 
angles at q — (ferromagnetic) and the step-like dashed line 
corresponds to q — 2k f (paramagnetic). For increasing r a the 
optimal spin-spiral wave vector becomes smaller, such that 
the Fermi spheres, separated at q — 2k f, begin to overlap. In 
order to gain energy the spin-up and spin-down orbitals in the 
overlapping region hybridize and the orbital angle 8 describes 
the mixing of the spin-up and spin-down states. 



(cf. Fig. [5]) also for the unoccupied states. In a comple- 
mentary work we have investigated the SSDW state using 
the Optimized-Effective-Potential (OEP) method within 
the framework of non-collinear Spin-Density-Functional 
Theory (SDFT)i&. In contrast to our findings within the 
OEP-DFT framework, i.e. , an effective single-particle 
theory restricted to local external potentials, here we do 
not find holes below the Fermi surface (cf . Ref . [l6| for de- 
tails). This is expected because it was shown in Ref. [l3 
that the HF ground state has no holes below the Fermi 
surface if the interaction is repulsive. Therefore our as- 
sumption of occupying only one band is justified. 

At the single-particle level we have an intuitive un- 
derstanding of the instability: as the two distinct spin- 
up and spin-down regions of the paramagnetic state are 
squeezed into each other, the orbitals in the overlapping 
region hybridize. This hybridization then leads to the 
opening of a direct gap between the HF single-particle 
dispersions corresponding to b — 1, 2 at k z = as well as 
to a lowering of both the symmetry and the total energy 
of the system. The mixing of the spin-up and spin-down 
orbitals is given by the orbital angles 8k, capable of de- 
scribing a continuous transition between the paramag- 
netic and the ferromagnetic state [Eqs. (|25a[) and (|25b[) 
respectively] . The behavior of the orbital angles at the 
optimal spin-spiral wave vector is shown in Fig. 



q/k f 

FIG. 8: (Color online) Total energies per electron of the 
SSDW state described with the density-matrix-power func- 
tional as a function of the spin-spiral wave vector q for vari- 
ous values of a at r 3 — 5.0. The total energy per electron at 
q = is subtracted in order to emphasize the behavior with 
increasing q. 



B. Correlated Functionals 

The density-matrix-power functional reduces to the 
uncorrelated HF approximation for a — 1 and to the 
Miillcr functional for a — 0.5. The latter one is known-^ 
to over-correlate and therefore one expects that decreas- 
ing a from 1 — * 0.5 increases the amount of correla- 
tion in the system. This picture was verified in Ref. [^, 
where an optimal value of a w 0.6 was found in the re- 
gions of metallic densities for the paramagnetic UEG. In 
Fig. [H the dependence of the total energy per particle 
at r s — 5.0 is shown for various a. It should be noted 
that the configuration for q > 2k f cannot be interpreted 
as the paramagnetic state in the correlated case. This 
is due to the fact that correlations smear out the sharp 
step found for the uncorrelated case in the momentum 
distribution around the Fermi surface (see Ref. 15 for de- 
tails). Therefore at q = 2k f the (fractionally) occupied 
regions in momentum space are not necessarily disjoint. 
Only when the occupied regions separate into two parts 
the configuration may corresponds to the paramagnetic 
state. However, the configuration at q — may still be 
interpreted as the ferromagnetic state. 

From Fig. [5] it is clear that the instability w.r.t. a 
SSDW is still present for a = 0.9. For higher values of a 
the instability disappears and for a = 0.5, 0.6 the energy 
has a maximum in the SSDW region. Thus for values of 
a which provide good correlation energies for the UEG 
in the paramagnetic regime there is no SSDW formation. 
In order to understand the reason for this it is instruc- 
tive to look at various contributions to the total energy. 
In Fig. [5] we compare the correlation energy contribution 
with the contribution coming from the kinetic and ex- 
change terms. The minimum is still present considering 
only kinetic and exchange contributions, but for decreas- 
ing a the correlation contribution damps out the insta- 
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FIG. 9: (Color online) Comparison of the correlation energy 
with the contribution from kinetic plus exchange terms (KEX) 
for a = 0.6, 0.7, 0.8 at r s = 5.0. All energy contributions are 
shifted such that the value at q = is zero. For decreasing a 
the minimum in the KEX contribution shifts to higher values 
of q. The correlation contribution however damps out this 
instability for values of a that yield good total energies at 
metallic densities. 



tions we define the relative correlation energy Se c as 



8e c r 



(30) 



In Fig. [TU] we show the dependence of this quantity on 
the spin-spiral wave vector for the correlation parame- 
ter a = 0.6. The absolute value of the relative correla- 
tion is smallest in the region of the SSDW instability 
(q = kf — > 2k f), which explains why the instability is no 
longer present when correlations are included. Further- 
more we can see that the relative correlation is dom- 
inant in the region of the ferromagnetic configuration. 
This can be understood by noticing that the density- 
matrix-power functional approximates the correlation en- 
ergy by a prefactor times a Fock integral (most present- 
day functionals in RDMFT approximate correlations in 
this wayi^^^^^^^*^). Since Fock integrals 
imply that equal spins are particularly correlated, one 
would expect a similar dependence of the relative corre- 
lation energy for other RDMFT functionals. 




i i i i i i i i i i i 

0.5 1 1.5 2 2.5 

q/k f 



FIG. 10: (Color online) Relative correlation energy as a func- 
tion of the spin-spiral wave vector q for various r s at a — 0.6. 
Correlations are more important around the ferromagnetic 
configuration. In the region of the Hartree-Fock SSDW insta- 
bility correlations have the smallest effect. This indicates why 
the instability is not sustained when correlations are included 
at the level of the density-matrix-power functional. 



bility more and more. One might suspect that at high 
densities, where exchange dominates correlations, the in- 
stability sustains. Our findings in Sec. IIV Al show that 
in the HF approximation the energy gain decreases when 
the density increases, which is consistent with an ana- 
lytic argument that at high densities the energy gain 
by forming a SDW and/or CDW is overcome by correla- 
tions. Furthermore our results indicate that correlation 
effects dominate the SSDW instability also at intermedi- 
ate densities. 

In order to gain some insight into the role of correla- 



V. SUMMARY AND CONCLUSION 

We have investigated the instability of the uniform 
electron gas w.r.t. the formation of a spin-spiral density 
wave within Rcduced-Density-Matrix- Functional Theory, 
which includes the Hartree-Fock approximation as an im- 
portant limiting case. To our knowledge this is the first 
numerical Hartree-Fock study of the spin-spiral state in 
the electron gas, despite the fact that Overhauser pre- 
sented his analytical work on the problem more than 
four decades ago. In Overhauser's work, the optimal 
spin-spiral wave vector was not determined. Our study 
shows that, in contrast to common belief, the optimal 
spin-spiral wave vector is not always close to 2k f. While 
at high densities we confirm this value for the optimal 
wave vector, for lower densities (just before the transi- 
tion to the ferromagnetic state) the optimal wave vectors 
even approaches kf. 

Within the framework of Reduced-Density-Matrix- 
Functional Theory we also studied the effect of corre- 
lations on the spin-spiral density wave instability using 
the recently proposed density-matrix-power functional. 
Not unexpectedly, we find that the inclusion of corre- 
lations suppresses the instability, which is explained by 
the behavior of the correlation energy in the region of the 
spin-spiral density wave instability. 
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